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Abstract. Finite volume methods for problems involving second order operators with full 
diffusion matrix can be used thanks to the definition of a discrete gradient for piecewise constant 
functions on unstructured meshes satisfying an orthogonality condition. This discrete gradient 
is shown to satisfy a strong convergence property on the interpolation of regular functions, and 
a weak one on functions bounded for a discrete H 1 norm. To highlight the importance of both 
properties, the convergence of the finite volume scheme on a homogeneous Dirichlet problem 
with full diffusion matrix is proven, and an error estimate is provided. Numerical tests show the 
actual accuracy of the method. 
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1 Introduction 

The approximation of convection diffusion problems in anisotropic media is an important issue in 
several engineering fields. Let us briefly review four particular situations where the discretization 
of a nondiagonal second order operator is required: 

1. In the case of a contaminant transported by a one-phase flow, one must account for the 
diffusion-dispersion operator div(AVtt), where the matrix A(x) = A(x)bj + ^(x)q(x) • q(x)* 
depends on the space variable x and q(x) is the velocity of the fluid flow in the porous 
medium. The real parameter \{x) corresponds to a resulting isotropic diffusion term, 
including dispersion in the directions orthogonal to the flow, and the real parameter fj,(x) 
to an additional diffusion in the direction of the flow The term q(x) is then given 
by q(x) = K(x)Vp(x), where p(x) is a pressure and K(x) another nondiagonal matrix 
(the absolute permeability matrix, depending on the geological layers), and satisfies the 
incompressibility equation divq(x) = 0. In this coupled problem, one must simultaneously 
compute this pressure and the contaminant concentration u(x). 

2. In the study of undersaturated flows in porous media (for example, air- water flows), two 
equations of conservation have to be solved, associated with two unknowns, pressure and 
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saturation. These equations include nonlinear hyperbolic and degenerate parabolic terms 
with respect to the saturation unknown. As in the preceding case, one must discretize such 
terms as divq(x) = div(K(x)Vp(x)), where again K(x) is a nondiagonal matrix depending 
on the geological layers. 

3. In the case of the compressible Navier-Stokes equations, one has to discretize the viscous 
forces operator, which can be written under the form aAu + 6Vdivu (a and b are deduced 
from the dynamic viscosity coefficients and u is the fluid velocity). In this problem, the 
term Vdivu involves all the cross derivatives <9? u. 

4. Some problems arising in financial mathematics lead to anisotropic diffusion equations 
in high-dimensional domains (dimension equal to 5 or more for example). Under some 
assumptions on financial markets [23]) the price of a European or an American option is 
obtained by solving a linear or nonlinear partial differential equation, involving the second 
order anisotropic diffusion matrix A = EE*, where E is a real matrix. 

All these cases involve a term under the form div(AVu), where A is a (generally) nondiagonal 
matrix depending on the space variable and u is a function of the space variable in steady 
problems, and of the space and time variables in transient problems. Finite element schemes are 
known to allow for an easy discretization of such a term on triangular or tetrahedral meshes |22j . 
However, in engineering situations such as the ones described above, one also has to discretize 
convection and reaction terms, and avoid numerical instabilities. Unfortunately, finite element 
methods (and more generally centered schemes) are known to generate instabilities on coarse 
grids, although some cures may be proposed, see |14l l3*]: therefore a great many numerical codes 
P] 01 HHEUE2] use finite volume or finite volume - finite element type schemes, which allow the 
implementation of discretization techniques (such as the classical upwind schemes) which prevent 
the apparition of instabilities. Let us also note that finite volume schemes are known for their 
simplicity of implementation, particularly so when discretizing coupled systems of equations of 
various nature. 

Besides, a thorough mathematical analysis has now been improved, showing that finite volume 
methods are well suited and convergent for a simple convection diffusion equation in the case 
where A(x) = X(x) Id- Indeed, this analysis has been completed (see |17j . |24j . jlHj . [H] ) in the 
case of grids (called admissible in the sense of [S] , see also Definition 12.11 below) satisfying an 
orthogonality condition: the line joining two cell centers is orthogonal to the interface between 
the two cells, thus ensuring a consistency property when approximating the normal flux at the 
cell interface by centered finite differences. Some examples of such admissible grids are the 
Delaunay triangular meshes or tetrahedral meshes, rectangular or parallelepipedic meshes in 2 
or 3 dimensions, and the Vorono'i meshes in any dimension. 

But the situation is quite different in the case where the condition A(x) = X(x) Li no longer 
holds: only few of the actual discretization methods used for handling nondiagonal second order 
terms on finite volume grids meet a full mathematical analysis of stability or convergence. Let 
us briefly review some of them. A first one, in the case where A(x) = A(x) M, where M is 
a symmetric positive definite matrix, consists in adapting the above orthogonality condition 
by stating that the line joining two cell centers is orthogonal to the interface between the two 
cells with respect to the dot product induced by the matrix A -1 . Indeed, it is also possible 
to consider the case where M depends on the discretization cell, by using, in each cell, the 
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orthogonal bisectors for the metric induced by M~ l (see ^Sl an d [H] section 11 page 815). In 
the case of triangular grids, this yields a well defined scheme under some restriction on the 
allowed anisotropy for a given geometry, since the cell center is chosen as the intersection of the 
orthogonal bisectors of the triangle for the metric defined by M~ l . Another method consists in 
defining the finite volume method as a dual method to a finite element one (for example, a PI 
finite element |2j or a Crouzeix-Raviart one, see e.g. |13j ) . 

Another possibility to derive a finite volume scheme on problems including anisotropic diffusion 
is to construct a local discrete gradient, allowing to get, at each edge a of the mesh, a consistent 
approximate value for the flux f (A(x)Vu(x)) ■r\. (T 6rf{x) involved in the finite volume scheme (n CT 
is a unit vector normal to the edge a, and d7(a;) is the d— 1 Lebesgue measure on the edge a ). In 
two space dimensions, such a scheme was introduced in <Jo on arbitrary meshes, but the proof of 
convergence was only possible on meshes close to parallelograms. Still in 2D, a technique using 
dual meshes is introduced in |191 [7j , which generalizes the idea of |25l I2(J| for div-curl problems 
to meshes with no orthogonality conditions; however the use of a dual mesh renders the scheme 
computationally expensive; moreover it does not seem to be easily extended to 3D. In |1U| . we 
used Raviart-Thomas shape functions, generalized to the case of any admissible mesh (again in 
the sense precised of (Hj, see also Definition 12 . 1 1 below) . in order to define a discrete gradient for 
piecewise constant functions. The strong convergence of this discrete gradient was then shown 
in the case of the elliptic equation —Au = f. A drawback of this definition was the difficulty 
to find an approximation of these generalized shape functions in other cases than triangles or 
rectangles. 

We therefore propose in this paper a new cheap and simple method of constructing a discrete 
gradient for a piecewise constant function, on arbitrary admissible meshes in any space dimension 
(this method has been first introduced in JJ). We prove that the discrete gradients of any 
sequence of piecewise constant functions converging to some u G Hq(Q) weakly converges to Vii 
in L 2 (J7). Moreover, the discrete gradient is shown to be consistent, in the sense that it satisfies 
a strong convergence property on the interpolation of regular function. In order to show the 
efficiency of this approximation method, we use this discrete gradient to design a scheme for the 
approximation of the weak solution u of the following diffusion problem with full anisotropic 
tensor: 

-div(AVu) = / in fi, , . 

u = on dU, ' 

under the following assumptions: 

£1 is an open bounded connected polygonal subset of M. d , d G N*, (2) 



A is a measurable function from O to Aid 
where A4d(M) denotes the set of d x d matrices, 

such that for a.e. A(x) is symmetric, , . 

and the set of its eigenvalues is included in [a(x), j3{x)\ 
where a,(3 £ L°°(Q) are such that 
< ao < a(x) < (3(x) for a.e. x G f2, 

and 

feL 2 (Q). (4) 
We give the classical weak formulation in the following definition. 
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i / A(x)Vu(x) ■ Vv(x)dx = [ f(x)v(x)dx, W € JZo(fi). 




Remark 1.1 For i/te sake of clarity, we restrict ourselves here to the numerical analysis of 
Problem f?)). however, the present analysis readily extends to convection- diffusion-reaction prob- 
lems and coupled problems. Indeed, we emphasize that proofs of convergence or error estimate 
can easily be adapted to such situations, since the discretization methods of all these terms are 
independent of one another, and the treatment of convection and reaction term is well-known 



The outline of this paper is the following. In Section |2 we present the method for approximating 
the gradient of a piecewise constant function, and we show some functional properties which help 
to understand why the present definition of a gradient is well suited for second order diffusion 
problems. In Section |31 we present the finite volume scheme for Problem ([I]), and we show the 
strong convergence of the discrete solution and of its discrete gradient. In Section @J we give 
an error estimate for Problem (fT|). and we illustrate this study by some numerical examples in 
Section |SJ Some short conclusions are drawn in Section HO 

2 A discrete gradient for piecewise constant functions 

We present in this section a method for the approximation of the gradient of piecewise constant 
functions, in the case of grids satisfying some orthogonality condition as defined below. 

2.1 Admissible discretization of f2 

We first present the following notion of admissible discretization, which is taken in jHj. The 
notations are summarized in Figure ^ for the particular case d = 2 (we recall that the case d > 3 
is considered as well). 



exact(see TWij or J$/). 




Figure 1: Notations for a control volume K in the case d = 2 
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In the following definition, we shall say that a bounded subset of M. d is polygonal if its boundary 
is included in the union of a finite number of hyperplanes. 

Definition 2.1 [Admissible discretization] Let O be an open bounded polygonal subset of 
R d , and 30 = \ O its boundary. An admissible finite volume discretization o/O, denoted by 
V, is given by V = (Ai,S,'P), where: 

• M is a finite family of non empty open polygonal convex disjoint subsets o/O (the "control 
volumes") such that = UxeM^. For any K G M., let dK = K\K be the boundary of 
K and m(K) > denote the measure of K . 

• £ is a finite family of disjoint subsets of O (the "edges" of the mesh), such that, for all 
a G £, there exists a hyperplane E of M, d and K G A4 with a = dK n E and a is a non 
empty open subset of E. We then denote by m a > the (d— 1)- dimensional measure of a. 
We assume that, for all K G M, there exists a subset £^ of £ such that dK = ^3o^E K ~v- It 
then results from the previous hypotheses that, for all a G £, either a C dQ or there exists 
(K, L) G M 2 with K L such that K n L = a; we denote in the latter case a = K\L. 

• V is a family of points of VL indexed by M., denoted by V = {xk)k^m- The coordinates 

(i) 

of xk are denoted by x K , i = l,...,d. The family V is such that, for all K G M, 
xk G K. Furthermore, for all a G £ such that there exists (K, L) 6 Ai 2 with a = K\L, it 
is assumed that the straight line (xk,%l) going through xk and xl is orthogonal to K\L. 
For all K G M and all a G £k, 1st z a be the orthogonal projection of xk on a. We suppose 
that z a G a if a C d£l. 

The following notations are used. The size of the discretization is defined by: 

hx> = sup{diam(_fC), K G A4}. 

For all K G M and a G £k, we denote by n^ jCr the unit vector normal to a outward to K. We 
denote by dx,a the Euclidean distance between xk and a. We then define 

T~K,c 



dK,a 

The set of interior (resp. boundary) edges is denoted by £ m t (resp. £ C xt), that is £ m t = {a G £; 
a <£_ 30} (resp. £ cxt = {c G £; a C 8ft}). For all K G M, we denote by Afx the subset of M 
of the neighbouring control volumes, and we denote by £K,ext = £k H £ ex t- For all a G £i n t> let 
K,L G M. be such that a = K\L; we define by d^\L the Euclidean distance between xk and 
xl, by x\kl the unit normal vector to K \L from K to L, and we set 

(6) 



For all a G £ ex t, let K £ M be such that a G £k\ w e define 

T a = T K ,a- (7) 

For all K G M. and a G £k, we define 

D K ,o = {tx K + (1 - t)y, t G (0, 1), y G a}, 
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For all a G £j n t, let K, L £ M be such that a = K\L; we set = -D_r- >0 - U Z?l,o-- For all cr G £ cx t , 
let K G M. be such that a G £k', we define D a = Dk,ct- 
For all a € £ , we define 

x CT = — p— r / x d7(x). (8) 
We shall measure the regularity of the mesh through the function 6t> defined by 



I'D 



(9) 



Definition 2.2 Let f2 fre an open bounded polygonal subset of M. d , and T> an admissible dis- 
cretization oftl in the sense of Definition \2.1\) . We define Hx> as the set of functions u G L 2 (0) 
which are constant in each control volume. For u G Hd, we denote by uk the constant value of 
u in K. We define the interpolation operator P-p ■ C(f2) — > Hx>, by u ^ Pt>u such that 

Pvu{x) = u(x K ) for a.e. x G K, \/K G M. (10) 

For (u,v) G {Hv) 2 an d f or an V function a G L°°(f2), we introduce the following symmetric 
bilinear form: 

[it, V ]x>,a = 

K\L&e int KeM o-G^A'.oxt 



where we set 

m _ = 



7— r / a(x)dx, Vcj G £. (12) 



Remark 2.1 One could also take, for a a , the harmonic averaging of the values in K and L 
when a = K\L. 

We then define a norm in Hx> (thanks to the discrete Poincare inequality (|13jl given below) by 

||w||x> = ([«, u] Vt i) 1/2 

(where 1 denotes the constant function equal to 1). Indeed, the discrete Poincare inequality 
writes (see [S]): 

IMIl 2 ^) ^ diam(fi)||u;||x>, Vw G Hd- (13) 

Let us now give a relative compactness result, which is also partly stated in some other papers 
concerning finite volume methods [S], |12) . 

Lemma 2.1 (Relative compactness in L 2 (Q)) Let f2 be an open bounded connected polygo- 
nal subset o/R , d G N* and let (T> n ,u n ) ne ^ be a sequence such that, for all n G N, T> n is an 
admissible finite volume discretization of Q in the sense of Definition \2.f\ and u„ G Hj) n {Q) (cf 
Definition \2.2\) . Let us assume that linij^oo hx>„ = 0, and that there exists C\ > such that 
\\u n \\v n < CfTJ, for all n G N. 

Then there exists a subsequence of (T> n ,u n ) n <^, again denoted (V n ,u n ) n ^, andu G Hq(Q) such 
that u n tends to u in L 2 (Q) as n — > +oo, and the inequality 

/ |Vn(x)| 2 drE < liminf \\u n \\ v (14) 
Jn n ^°° 
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holds. Moreover, for all function a G L°°(0), we have 



lim [un, PT> n <p]v n , 



a(x)Vu(x) ■ Vip(x)dx, Vtp G C~(fi). 



(15) 



Proof. The proof of the existence of the subsequence again denoted (V n ,u n ) ne ^, and of 
u G Hq(Q) such that u n tends to u in L 2 (J7) as n — > oo, is given in [S]. Assertion (|14[) was 
proven in ^21 (Lemma 5.2). Let us first show 1)15(1 in the case a G C 1 (J7). Let G C£°(f2). 



Defining, for all n G N, Tf 
lim Bfr) 



(«) 



f^u n (x)div(a(x)'Vip(x))dx, we get that 



u(x)div(a(x)V tp(x))dx 



a(x)Vu(x) ■ Vip(x)dx. 



We consider a value n sufficiently large such that for all K G M n and x G if , if 7^ then 



Mr n 9^ = 0. Defining T 2 (n) = [u n , PvMv n 



Tj^' we obtain 



(»0 



El 



XI m(if|L)(uL - u k )Rkl, 



with 



a(x)Vy>(a;) • n KL dj{x), VK eM, VL G Af K - 



K\L 



Since there exists some real value C2, which does not depend on T> n , such that < ^f 1 ©,^ 

we conclude in a similar way as in jHj that lirrin^oo = 0, which gives (|15|) in this case. Let 
us now consider the general case a G L°°(0,). Let e > be given. We first choose a function 
a G C 1 (0) such that \\a — ^ £ - Then we have, for all n G N, using the Cauchy-Schwarz 

inequality, 

([U n , PT>Mv n ,a ~ [Un, pD n <Ph n ,a) 2 < ^ T K \ L (a K L ~ <*Kl) 2 I <p(x L ) ~ <p(x K )\ 2 

K\Le£ int 

x ^2 t~k\l\ u l ~ u K \ 2 

K\Le£ int 

and therefore, setting C3 = ||Vv?||.L<»(n)j the properties \ip(x£)— <p{xk)\ < C^k\l &ndm(K\L)d K \ L 
d m(D K \ L ) lead to 



(K, P-D n <p]v n ,& - [u n , PD n <p]v n ,a) < d Qffl 
In the same manner, we get 



h *a-a\\l2 {n) Ofa< d qgp 2 qrj 



a(x)Vu(x) ■ Vtp(x)dx 



J a(x)Vu(x) ■ V<p(x)dx\ < C^ 2 \\Vu\\l w . 



Since a G C 1 (0), we can apply (|15|). proven above for such a function. It then suffices to choose 
n large enough such that 



to prove that 



[Un, Pv n <pfo n , 



a(x)Vu(x) ■ Vif{x)dx < e, 



a(x)Vu(x) ■ V(p(x)dx < C4 e, 



where the real > does not depend on n. This concludes the proof of (|15|) in the general 
case. □ 



7 



2.2 Definition of a discrete gradient 

We now define a discrete gradient for piecewise constant functions on an admissible discretiza- 
tion. 

Definition 2.3 (Discrete gradient) Let be an open bounded connected polygonal subset of 
W 1 , d £ N*. Let T> = (A4,£,V) be an admissible finite volume discretization of VL in the sense 
of Definition ^. 11 Let us define, for all K £ A4, for all L £ Mk, 

A k ,l = t k \ l (x k \ l - x K ), (16) 

and for all a 6 £K,exk> we define 

A K ,a = r a (x a - x K )- (17) 
We define the discrete gradient Vx> : H?> — ► H^, for any u £ Hx>, by: 
Vt>u(x) = (Vt>u) k 

= Y Ak > l ( u l-u k )- Y Ak >° Uk 

y ' \LeAf K o-e£jr,ext 
for a.e. x £ K, MK £ M. 

Let us first state a bound for the L 2 (Q) d norm of the discrete gradient of any element of ifx>. 

Lemma 2.2 (Bound for Vptt) Let Q be an open bounded connected polygonal subset of M. d , 
d £ N* ; let T> be an admissible finite volume discretization of £1 in the sense of Definition \2.f\ 
and let 6 £ (0, 6-p]- Then, there exists C$, only depending on d and 6, such that, for all u £ Hx>: 

\\Vvu\\ L 2 {n) d < q^Mix). (is) 

Proof. Let u £ Hjy- Let us denote, for all K £ Ml, L £ Mk and a = K\L, Sk^u = ul — uk, 
and for a £ £K,ext, b~K,aU = —uk- Then Definition (TTTj) leads to 



A 

and Definition (|2.3|) leads, for a given K £ A4, to 



Y I \ Y T K\l($K,K\Lu) 2 + Y T °( 6 K,aU) 2 J , 



Vn{K){V V u) K = Y T <r( x a - x k)8k,<tU. 
a&£ K 



Using the Cauchy-Scharwz inequality, we obtain 

m^KV^I 2 ^ Y, I 2 Y T a( S K,aU) 2 , 

and, since, for a £ £k, one has \x a — xk\ = d(x<7,xx) < "f 12 "! 

m(K) 2 \(V v u) K \ 2 < ^H°)d K ,o Y T -( 6 K,au) 2 . (19) 



8 



Since YlaeE K M a ) d K,a = d m(K), (JHJ) gives: 

d 



e 2 

a££ K 



Summing over K G M, we get 

||V7?«||| a(n)< , < ^\\u\\l- 
which gives (JTHJ with (M = ■ □ 

We now state a weak convergence property for the discrete gradient. 
Lemma 2.3 (Weak convergence of the discrete gradient) 

Let Q, be an open bounded connected polygonal subset ofM. d , d G N* ; let T> be an admissible finite 
volume discretization of £1 in the sense of Definition \2. A We assume that there exist ux> G Hx> 
and a function u € Hq(Q) such that u?> tends to u in L 2 (J7) as hx> tends to while \\ud\\t> 
remains bounded. Then V^hd weakly tends to V« in L 2 (Q) d as hx> — ► 0. 

Proof. Let ip G C£°(0). We assume that hx> is small enough to ensure that for all K G M 
and x G K, if ip(x) / then £x,ext = 0- The expression T^, defined by 

= I Pv t p{x)Vx)U T ,(x)dx, 

satisfies, using (JTJJJ), 

2j§ = ^2 t k\l(ul ~ u K ) {{x k \l ~ x K )ip(x K ) + (xl ~ x K \ L )ip(x L )) , 

K\Le£ int 



where we denote, for the sake of simplicity, uk = (ut>)k f° r all X G Ai. We thus get 
If + Tf with 

2g| = 2^ t k\l{ul-u k ){xl-x k ) 

and 

2j§ = X] t k\l(ul-uk){x k \l 2 - ^O 1 *))- 

Thanks to the Cauchy-Schwarz inequality, we get 

( :Z fi) 2 - Yl t k\l{ul-u k ) 2 t k\l{w{xl) - ip{x K )) 2 \x K \ L _ XL+ ^ XK |2. 

Since — < ^l^xli ~~ x il + \\ x k\l ~ x k\ < /tx>, there exists 6*6 > 0, depending on 

d, and ip such that, 

(^g) 2 < \\ U v\\lqg%m(n), 

and therefore we get 



lim 2^ = 0. 
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We then compare 2gj with 

Tf = - / u v (x)V(p(x)dx = (ul~ u K ) \ (p(x)n KiL d'j(x). 

JO. J K\L 



Since 



and since 



n K ,L 



K\Lee in 

XL - X K 



if(x)d , y(x) 



dx\L 
<p(x K ) + tp(x L ) 



HK\D Jk\l 
we get, thanks to the Cauchy-Schwarz inequality, 

Since 



< \\Vtp\\ L °°(n) hv, 



u(x)Vip(x)dx = / tp{x)V u{x)dx , 



lim T§ = - 

we have thus proven, thanks to the density of C£°(fi) in L 2 (ft), the weak convergence of Vt>ut> 
to V«(x) as hx> 0. This completes the proof of the lemma. □ 

We now study, for a regular function ip, the strong convergence of the discrete gradient Vx>Pvf 
to V93. This study uses the following lemma. 

Lemma 2.4 Let f2 be an open bounded connected polygonal subset ofM. d , d G N*, let T> be an 
admissible finite volume discretization of Q in the sense of Deftnition \2.1\ Then we have 



1 



m(K) 



m(a)(x a - xo) (iiK,a- • «), VK G M, Vx G R d , Vu G 



(20) 



Proof. For any K G A4, we denote, for a.e. x G <9-KT, by ng^(x) the normal vector to dK 
at the point x outward K. Let v and w £ R d be given. We have, considering vectors as d x 1 
matrices, and denoting by w t the transposed 1 x d matrix of w, 



w 



dK 



(x - x )n t K (x)d'y(x) ) v 
dK J 

uf (x — xq) v t riK(x)dj(x) = 



w l (x — xq) n K (x) v dj(x) 
dK 

(v (x — xq) 1 w) ■ nx(x)d'y(x) 



dK 



div(v (x — xq) 1 w)dx = m(K) v l w. 



K 



This gives fRft. □ 



Lemma 2.5 (Consistency property of the discrete gradient) Let £1 be an open bounded 
connected polygonal subset ofM. d , d G N*, let V be an admissible finite volume discretization in 
the sense of Definition \2.l\ and let 9 G (0, 9x>]- Let u G C 2 (f2) be such that u = on the boundary 
ofVt. Then, there exists Cj, only depending on Vt, 9 and u, such that: 



\\V v Pvu - Vu\\ L a(ny < h v . 
(Recall that Px> is defined by 110]) and Vp in Definition \2.'A ) 



(21) 
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Proof. From Definition 12.31 and 1)1(1 1) . we can write for any K G Ai 
m.{K)(V v P v u) K = E t k\l(xk\l ~ x K )(u{x L ) - u{x K )) - E T a (x a - x K )u(x K ). (22) 

Let (¥u)k be the mean value of Vu on K: 

(Vu) K = / Vu(x)dx. 

Thanks to the regularity of u (and the fact that u = on the boundary of $1), there exists C%, 
only depending on u (indeed, Cjgjonly depends on the L°°-norm of the second derivatives of u), 
such that, for all a = K\L G £i n t, 

w(x^) — U\X ft) 

|e CT | < Cj^fiD, with e a = {Vu) K ■ -j , (23) 



and, for all a E £kt 



ext j 



|e CT | < Cjgfto, with e a = (Vu) K • n^ j(J . (24) 



Thanks to (|22 |) . and ijM jl . we get, for all K £ M: 

m(K)(VvPvu)K = E m (°")(^ - ^)(Vu)k • n Kj(J . + ifo, 

with Rk = — e cr m ( (7 ) < ^( x cTi x K)- Applying (|2"U)l gives 

<t<=£k 

m(K)(V v P v u) K = m(K)(Vu) K + Rk- (25) 
Using the inequalities (|2"5")l and (|24[). we have 

|^| < E m ^)^- = ^^m(K). (26) 

Then, from ()25() and (|26[). we obtain 

E |(VpP c ?% - (Vnk|V(iO < 

/^ffe,\ 2 „ ( 2 7) 



E (^) 4- 



In order to conclude, we remark that, thanks to the regularity of u, there exists Cg, only 
depending on u (here also, only depends on the L°°-norm of the second derivatives of u), 



such that: 



E f \Vu(x)-(Vu) K \ 2 dx<G^hl. (28) 



Then, using l)27|) and (|28|). we get the existence of Cjj|, only depending on f2, and u, such that 
J2U) holds. □ 
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Remark 2.2 (Choice of the points xk and x a ) Note that in the vroof of Lemma \2.c\ one is 

free to choose any point lying on K\L instead of x^\l * n the definition of the coefficients Ak,l- 
However, we need this choice in the proof of the strong consistency of the discrete gradient 
(Lemma \2. 5\) . Conversely, in the proof of Lemma \2. 51 we could take any point of K instead of 
xk in the definition of Ax^l- However, the choice of xx is crucial in the proof of Lemma \2.,'Jl 
when comparing the terms 2pa and 1fc\ one needs the property of consistency of the normal flux, 



which follows from the fact that t&k l ~ 



Lemma 2.6 (A sufficient condition for the strong convergence of the discrete gra- 
dient) 

Let Q be an open bounded connected polygonal subset ofM. d , d G N*, let 6 > and let T> be an 
admissible finite volume discretizations in the sense of Definition \2.l[ such that 9x> > 0. Assume 
that there exists a function u-p G Hx> and a function u G Hq{Q) such that ux> tends to u in 
L 2 (J7) as hp tend to 0. Assume also that there exists a function a G L°°(J7) and Qo > such 
that a(x) > ao for a.e. x G Q and [ux>, u-p]v,a tends to J^a(x)Vu(x) 2 dx as hx> tends to 0. 
Then Vx>uv tends to Vu in L 2 {Q) d as hx> tends to 0. 

Proof. Let ip G C£° (f2) be given (this function is devoted to approximate u in Hq{Q)). Thanks 
to the Cauchy-Schwarz inequality, we have 



with 



and 



[ (V v uv(x) - Vu(x)) 2 dx < 3 (7f + Tf + T 9 ) 
Jn 

2§ = / (Vuud(x) " VvPv^(x)) 2 dx, 
T i = / (VvPv^(x) - V<p(x)) 2 dx, 

1fa= J (V<p(x) - Vu(x)) 2 dx. 



We have, thanks to Lemma 12.51 
Thanks to Lemma 12.21 we have 



h lim o 2g = 0. (29) 



/ (V v v(x)) 2 dx < q&[v t v\D,i < -^[»,«Ka, G H V . 
Jn m ao 



We thus get, setting v = ux> — Pv^P in the above inequality, that 

2 



2s2 < ■^-([uv,v^]v,a - 2[uv,Pv<p]T>, a + [Pd<P, PD¥>]v,a) 



We have, applying twice Lemma l2~Tl that 

lim [ud, Px>(p]v,a = / a{x)V u{x) ■ V (p{x)dx (30) 
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and 

lim [Pp<f, Pvv]v,a = / a(x)V(p(xfdx. (31) 
Under the hypotheses of the lemma, we then get that 

_ SB 



limsup^ < -EL / a(x)(Vu(x) - V<p(x)) 2 dx. 

u_ ,n LLI n/.n /o 



We then get, gathering the above results, setting C\q = ess sup^gQ cv(x) + 1, that 

f (V v uv(x) - Vu(x)) 2 dx < qini / (V^(x) - Vu(x)) 2 dx + t\ 
Jn — Jn 
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with 



lim Egg =0. (32) 



Let e > 0. We can choose ip such that / r2 (Vy?(x) — Vu(x)) 2 dx < e, and we can then choose hx> 
such that < e. This completes the proof that 



lim / (Vt)Ut>(x) — Vn(x)) 2 dx = 0. 



(33) 



□ 



Remark 2.3 Thanks to Lemma \2. 61 we get the strong convergence of the discrete gradient in the 
case of the classical finite volume scheme for an isotropic problem. Note that in the above proof, 
we did not use the weak convergence of the discrete gradient, and therefore any point of K can 
be taken instead of xk in the definition of the coefficients Ak,l ■ We thus find that the average 
value in K of the gradient defined in 1101/ is also strongly convergent (the average of this gradient, 
defined by the generalized Raviart- Thomas basis functions, is obtained by replacing xk by the 
barycenter of K in the definition of Ak,l)- Note that the drawback of the generalization of the 
Raviart- Thomas basis was the difficulty for computing approximate values of the gradients. This 
drawback no longer exists for an averaged gradient. Nevertheless, the properties of convergence 
of the finite volume method shown here for non isotropic problems are only proven for the choice 
111 6)) in the definition of Ak,l, and not for the Raviart- Thomas basis. 



3 Application to Problem (0Q) 

3.1 The finite volume scheme 

Under hypotheses let T> be an admissible discretization of O in the sense of Definition 

12.11 The finite volume approximation to Problem (^Q) is given as the solution of the following 
equation: 

f n {A(x) - a(x)l d )V v u v (x) ■ Vx>v(x)dx + [u v , v]x>, a = J n f(x)v(x)dx, Vu G H v , 

denoting by 1^ the identity application of M. d . The existence and the uniqueness of the solution 
ud to will be stated in Lemma 13.11 Note that in this formulation, we use the discrete 
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gradient on part of the the operator only, while on a homogeneous part, we write the usual cell 
centered scheme. This needs to be done in order to obtain the stability of the scheme, that is 
some a priori estimate on the discrete solution. If we take a = in (|34[). we are no longer able to 
prove the discrete H 1 estimate (j39j) below. Taking for v the characteristic function of a control 
volume K in (|34|). we may note that Equation (|34|) is equivalent to finding the values (uk)kgm 
(we again denote uk instead of {ud)k)i solution of the following system of equations: 



where 



2 Fkl+ Fk ° = [j( x ) dx i VKgM, 



(35) 



Fkl = T K \ L a K \ L (u K -u L ) + { A L A LK ■ Vx>u L - A K A KL ■ Vt>u k ) VK\L G & mt , (36) 

and 

Frg = T K aa a u K + A K A Ka ■ Vvur Vct G £x,cxt- (37) 
In (15*61 and (|37j). the matrices (A^)^ e _A^ are defined by: 

A ^ = — TIFT / ( A ( x ) " «(a;)Id)dx. (38) 
m{K) J K 

On can then complete the discrete expressions of Fkl and Fko using Definition 12.31 for Akl 
Akg-, and Vv^K for all if G .M, L G A/"^ and <r G f^. 
This is indeed a finite volume scheme, since 

Fkl = -F LK , Vif |L g £ int . 
The existence of a solution to will be proven below. 

3.2 Discrete H 1 ^) estimate 
We now prove the following estimate: 

Lemma 3.1 [Discrete H 1 estimate] Under hypotheses JU)-([""p, let T> be an admissible dis- 
cretization of £1 in the sense of Deflnition \2.1\ Let u G Hx> be a solution to \34\) - Then the 
following inequalities hold: 

ao\W\\v < diam(0)||/|| (X 2(o))2, (39) 
PROOF. We apply (jHU) setting v = u . We get 



(A(x) — a(x)ld)Vvu(x) ■ Vx>u{x)dx + [u, u\x> )a = / f(x)u(x)dx, 
n Jn 

which implies 

ao[u, u]-p < / f(x)u(x)dx. 
Jn 

Then the conclusion follows from the discrete Poincare inequality H13|) . □ 

We can now state the existence and the uniqueness of a discrete solution to (J3U). 
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Corollary 3.1 [Existence and uniqueness of a solution to the finite volume scheme] 

Under hypotheses let T> be an admissible discretization of Q in the sense of Definition 

\2.1\ Then there exists a unique u-p solution to \34\l - 

Proof. System (|34|) is a linear system. Assume that / = 0. Prom the discrete Poincare 
inequality (|13[). we get that u = 0. This proves that the linear system (|34|) is invertible. □ 

3.3 Convergence 

We have the following result, which states the convergence of the scheme (|34j) . 

Theorem 3.1 [Convergence of the finite volume scheme] Under hypotheses H^-jOU, let 
9 > 0. Let D be an admissible discretization of Q in the sense of Definition WIS, such that 
9x> > 9. Let uj) E Hd(Q) be the solution to \34\)- Then 



• ux> converges in L 2 (Q) to u, weak solution of Problem (GJ) in the sense of Definition M.A 

• the discrete gradient Vx>uv converges in L 2 (VL) d to Vu, 
as hjy tends to 0. 

Proof. We consider a sequence of admissible discretizations (£> n )neN such that hx> n tend to 
as n — > oo and 9x> n > 9 for all n £ N. Thanks to Lemma [3 .11 we can apply the compactness result 
(|2.1() . which gives the existence of a subsequence (again denoted (2? n ) ne pj), and of u E Hq{Q) 
such that ux> n (given by (|34() with T> = T> n ) tends to u in L 2 (0) as n —> oo. Let ip E C£° (O) be 
given, we choose v = P-v^ as test function in (|34|) . We obtain 

/ (A(x)-a(x)I d )Vi ?n ^ n (a;)-V2 5n P I?n ^(x)dx + [n I)n ,P I?n ^K )a = / /(x)P Dn ^(a;)dx. (40) 
Jn Jn 

We let n — > oo in 1)40(1 . Thanks to Lemma 12.31 and Lemma 12.51 (which provide a weak/strong 
convergence result), we get that 

lim / (A(x) - a{x)l d )V v n u Vn { x ) ' Vx^Pp^aOdx = / (A(x) - a(x)I rf )V-u(x) • V(/?(x)dx. 
n ->°°Jn Jn 

Using Lemma 12.11 we get that 

lim [uv n , Pv n <p]v n ,a = / a(x)Vu(x) ■ Vip(x)dx. 
n ^°° Jn 

Since it is easy to see that 

lim / f(x)Px> n (p(x)dx = / f(x)ip(x)dx, 
n ^°°Jn Jn 

we thus get that any limit u of a subsequence of solutions satisfies © with v = ip. A classical 
density argument and the uniqueness of the solution to © permit to conclude to the convergence 
in L 2 (Q) of ux> to u, weak solution of the problem in the sense of Definition ll.il as hx> tends to 
0, thanks to the fact that 9jy > 9. Let us now prove the strong convergence of Vpn© to VS. 
We have, using (|3"4"|) with v = ud, 

(A(x) - a(x)l d )V vu-d{x) ■ Vvuv{x)dx = / f(x)uv{x)dx - [u-d, «x>]x>,a- (41) 

Jn 



15 



Thanks to Lemma 12. 11 we have 

a(x)Vu(x) 2 dx < liminffu©, ud\v a, 



and therefore, passing to the limit in (|4*T|) . we get that 



limsup / (A(x) — a(x)I(i)Vpup(x) • Vj)Uv(x)dx < / f(x)ux>(x)dx — / a(x)Vu(x) 2 dx. 
Ht>— >-o Jn Jn Jn 



I'D 

We then have, letting v = u in © 



(A(x) - Q(x)I d )Vn(x) • Vn(x)dx = / f(x)u(x)dx - I a(x)Vu(x) z dx. (42) 
n Jn Jn 



This leads to 



limsup / {K{x) — ^x)!^)^ i)Uv{x) ■ V vuv{x)dx < I (A(x) — a(x)I^)Vn(x) • Vu(x)dx. 
h-v^o Jn Jn 

Using Lemma 12.31 which states the weak convergence of the gradient V-puv to Wu, we get that 

/ (A(x) — a(x)ld)Vu(x) ■ Vu(x)dx < liminf / (A(x) — a(x)ld)V vuv{x) ■ Vvuv{x)dx. 
Jn /i-d^o J n 

The above inequalities yield 



lim / (A(x) - a(x)I d )Vvuv(x) ■ Vvuv(x)dx = / (A(x) - a(x)l d )Vu(x) ■ Vu(x)dx. (43) 



From (JUJ), (j32) and (O, we thus obtain that 

lim [uv, ux>]v a = / a(x)Vu(x) 2 dx, 
h v^o Jn 



Therefore we can apply Lemma 12.61 This completes the proof of the strong convergence of the 
discrete gradient. □ 

4 Error estimate 

We now give an error estimate, assuming first that the solution of © is in C 2 (S7). In Theorem 
14.21 we will consider the weaker hypothesis that the solution of ((SJ) is only in H 2 (Q) under the 
assumption d < 3. 

Theorem 4.1 (C 2 error estimate) Assume hypotheses and that A and a are of class 

C 1 on $7. Let T> be an admissible finite volume discretization (in the sense of Definition \2.1\l . 
Let G (0, 9x>], where 6x> is defined by Let ux> G Hv be the solution of and u G Hq (£1) 
be the solution of |3J). We assume that u G C 2 (f2). 
Let us first assume that 

Va G £ cx t, / A(x)n dn {x) ■ (x a - z a )d 1 (x) = 0, (44) 

J a 
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where ngn(x) is the unit normal vector to d£l at point x, outward to fL 

Then, there exists C\\ only depending on £1, 9, olq, a, (3, A and ||?i||c 2 (Q); such that: 

\\u v - Pyu\\y < tfnfi-D, (45) 

\\uy - u\\ L 2 (n) < tfrftiy, (46) 

and 

\\Vyuy - Vu|| L 2 (n)d < Qrip-D- (47) 

Let us then assume that \44\ ) no longer holds, then there exists C\2, only depending onVt, 9, a, 
(3, A and \\u\\ H 2^, such that \7U\ ), J7i| ), j7ij| ) hold with C^J^/hy instead of Gffjpy- 

Remark 4.1 Let us give some sufficient (and practical) conditions for \44\l t° hold : 



• If the normal vector to dtt is an eigenvector ofA(x) for a.e. x G <90, then \44\) holds. Since 
this property is always satisfied in the isotropic case, the error estimate on the gradient 
j47\ ) holds for the classical cell centered scheme, for any admissible mesh. 

• If for all a G £ cxt with a G £k > the bary center x a of a is equal to the orthogonal projection 
z a of xk on a, then \44\) holds. This hypothesis is easy to ensure on rectangular and 
triangular meshes. 

Note also that one could replace \44\ l by \z a — x a \ < ^diam(K) (hy) 5 for all a G £ C xt- 

Proof. In the proof, we denote by Cj (i G N), various quantities only depending on $7, 9, a^, 
a, (3, A and ||u|| C 2( C ). 

Step 1. Let v G Hy. We first perform a computation of a consistency error, namely a bound 
for |Tn(v)| where Tffi^v) is defined by: 

/ (A(z) - a(x)l d )V v Pyu(x) ■ Vyv(x)dx + [Pyu, v]y, a = [ f(x)v(x)dx + atrfu). (48) 
Jn Jn 

We first consider the second term of the left hand side of (|48|) . Using classical consistency error 
(also used in the proof of Lemma 12,1(1 . one has: 

[Pvu,v]y, a = - dw{aVu)(x)v(x)dx + T 12 (v), (49) 
Jn 

with 

|2jt2fV)| < ^2m(a)\R a \6 a v, 

where 5 a v = \vk — vl\ if <t = K\L is an interior edge, 5 a v = \vk\ is a G £ cx t and \R a \ < C\%hy. 
Using the Cauchy-Schwarz inequality, this leads to: 

\^v)\ < C 14 hy\\v\\y. (50) 
We now consider the first term of the left hand side of (1481) . We have 



/ (A(x) - a(x)l d )VyPyu(x) ■ Vyv(x)dx = T 13 (v) + T 14 (v), (51) 
Jn 
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with 



Ifttfv) = J (A(x) - a(x)I d )Vu(x) ■ V v v(x)dx 



and 

IffillMI < C 15 \\Vt>Pvu - Vit|| i 2(Q)d||Vi5w||i2(n)<»- 
Using Lemma 12.51 and Lemma 12.21 we obtain 

tfEfv)\ < C w h v \\v\\ v . (52) 

We now compute ffig^)- For K £ M and a £ 8, let fix and respectively be the mean values 
of (A(a;) — a(x)Id)Vu on K and a: 

= 7T7x / (A(x) - a(x)Id)Vu(x)dz, /x CT = — !— ■ / (A(x) - a(x)I d Vu(x)d7(x). 

The regularity of u, A and a gives, for all K £ A4 and all a E (recall that | • | denotes the 
Euclidean norm in Mr): 

\HK ~ /v| < Cvjh-D- (53) 

Indeed, Q[T1 only depends on the L°°-norms of A, a and Vu and on the L°°-norms of the 
derivatives of A, a and Vtt. 

We now use ()53|) in order to give a bound of ^f^ v) as a function of /i£>. Indeed, the definition 
of Vx>v leads to: 

Tfn&v) = Y ^ K ' M k ){^vv)k = 
KeM 

v-K ■ a k>l (v L - v K ) - ^ a k>(T vk 

KeM \L&Af K o-e^A-.cxt 

X] Y ^K\L ■ A K ,L (VL ~ V K ) ~ Y A K,aV K ) +Ti 5 (v), 

KeM \LeAf K o-e£jf,ext 



with 



KeM \LeAf K o-G£x : oxt / 

qXT^ip Y (\ a k,l\ + \A l> k\)\vl-vk\+ Y l-^IK-l 

\a=K\Le£ int KeMcree K}C xt 

Since i/c,L = tk\l{ x k\l ~ x k) and Ax, a = 7"o-(£<t — one deduces from the preceding 

inequality, thanks to the definition of 9x> (which gives d(x a , xk) < {dx,a/&) if cr S £r) and using 
Cauchy-Schwarz Inequality: 

\^gv)\ < C 18 h v \\v\\ v . (54) 

We now remark that: 

IfT^v) ~ = Y Y ^ K \ L ' Ak,l (vl -vk)~ Y ^ ' A k,^ v k 

KeM \LeAf K <ree K}CXt J (55) 

Y ■ (x L ~ Xk)t c {v L - V K ) - Y Y ■ ( x * - X K )T a V K - 

a=K\Le£i n t KeM aeS Ktext 
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For a G £; nt , one has a = K\L and (xl — %k) = d a TiK,a where njf jCT is the normal vector to a 
exterior to K. 

For a G £ ex t) ° ne has °" £ £rt- Thanks to the fact that under homogeneous Dirichlet boundary 
conditions, the gradient of u is normal to the boundary, using Assumption (|44|) . we get that 



fi a ■ (x a - x K )T a = I (A(x) - a(x)l d Vu(x) ■ n dn (x)dj(x). 

J a 

Then, one deduces from l)55|) : 

#) " T^v) =-J div((A - al d )Vu){x)v{x)dx. (56) 

Therefore, since — div(AVn) = /, one has (|15|) with 2jjj](t>) = ffi^) + Tft^v) + ). This 
gives, with (jSn j) . (tBIl) : 

l^u)! < C7i 9 /H3||t;||2>. (57) 

This concludes Step 1. 
Step 2. 

Let ex> = Pt>u — ux> be the discrete discretization error. Using (|4"5)l and PH) give, for all v G -f/x>: 

y (A(x) - a(x)I d )V v e v (x) ■ V v v(x)dx + [e v ,v] v , a = ffilt^)- 

Taking u = ex> in this formula gives, with (|57j) . [e-D,ev]v,a < Qr9"f l x>ll e x , llx> and then, with 
C20 = QrW«o (since a ||e£>||x> < [ev,ev]v,a)- 

\\ev\\v < Qxfrv, (58) 

which is exactly ()45|) . 

Using the Discrete Poincare Estimate ()13|) and the fact that -u G C(O), one deduces (|46|l from 

(USD- 

The last estimate, Estimate (jlTj) . is a direct consequence of ((SHI), (pTj) and (|T%|) . This concludes 
the first part of the theorem, i.e. assuming (j44|) . 

If V no longer satisfies the hypothesis (f4"4"|) . one has to replace (|56|) by: 

- r ^jv) = ~ J div((A - aI d )Vn)(x)u(2;)dx + 2i 6 (u), 
where, recalling that by z a the orthogonal projection of xk on a (see Definition 12. 1(1 : 

= E E ■ ( z<j ~ ■ 

Thanks to the Cauchy-Schwarz inequality, we get 

#) 2 < E E T CT ^(diam(K)) 2 £ £ r CT 4, 
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which leads to 

< ^mcanjiivuiijjoiMiS,, 

i 

where m(d{l) is the d— 1-dimensional Lebesgue measure of dQ. This gives (|57)1 with instead 
of hx>. Following Step 2, this allows to conclude the proof. □ 

We now want an error estimate when the solution of ® is in H 2 (Q) instead of C 2 (f2), in the 
case where the space dimension is lower or equal to 3. Indeed, the C 2 -regularity of the solution 
of (JHJ) was used, in the preceding proofs, only four times, namely to prove (|2*3*|) . (|2*4"|) and (|2*H|) in 
Lemma 12.51 and to prove ()53j) in Theorem 14. II fin fact, it is also used for the classical consistency 
error (|49[). but, for this term, the generalization to the case where the solution of ((SJ) is in H 2 (Q) 
instead of C 2 (f2), in the case d < 3, is already done in We will now prove similar inequalities 
for u G H 2 (Q) n Hq(Q) (instead of u G C 2 (0) with u = on the boundary of Q) which will allow 
us to obtain the desired error estimate. 

Lemma 4.1 (Consistency of the gradient, u G H 2 (Q)) Under hypothesis with d < 3, 
let T> be an admissible finite volume discretization in the sense of Definition \2.1i and let 9 G 
(0,#x>]- Let u G H 2 (£l) n Hq(Q). Then, there exists C21, only depending on ft, 6 and u, such 
that: 

\\V v (P v u) - Vu\\ L 2 {n)d < (%T[pv\\u\\ H 2 {n) . (59) 
( Recall that Px> is defined in \10\i and Vx> in Definition \2.<A ) 

Proof. 

The proof follows the proof of Lemma 12.51 (in particular, recall that H 2 (Q) C C(f2) since d < 3). 
The C 2 -regularity was only used to prove (|23|) . (|24|) . (|28jl . We now prove similar inequalities in 
the case u G H 2 (Q). 

We begin with providing inequalities similar to (|23|) . ()24l) . We denote by (Vu) CT the mean value 
of Vu on a (recall that (Vu)k is the mean value of Vu on K). We use Inequality (9.63) of [H] 
(in the proof of Theorem 9.4, using the i^ 2 -regularity) . This inequality states the existence of 
C22, only depending on d and 9, such that, for all a = K\L G & in t: 

l^l 2 ^ G m^hr I l^)WI 2dz > with E - = (VH)„ • n Kl . - n(XL) 7 a(XK) , (60) 
and, for all cj G £ cx t ; if c G 

l^l 2 ^ °mrr%r I l# (^)l 2 dz, with e ct = (va) CT • n K , CT - (6 i) 

where: 

d 

|#(u)(z)| 2 = £ |AZ^(,)| 2 . 

We have now to compare (Vu) a and (Vw)^. This is possible thanks to Inequality (9.38) in 
Lemma 9.4 of [S]. Following this result, there exists C23, only depending on d and 9, such that, 
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for all K G M, all a £ £ K and all v G H 1 ^): 

f v(x)dx 3_ f v ( x )dj(x) < Cjtei diam ^ A - / \Vv(x)\ 2 dx < 



h 2 - r (62) 



Using (j62() with the derivatives of «, one deduces from 1)60(1 and l[61|). that there exists some real 
value C24 only depending on d and 6 such that 

^ qfea-rlr / l ff (*)(*)l 2 <k. with e - = ( v *)k • n K , CT - Q(xL) 7 a(xK) , (63) 

l ^ ±l m(cj)ii (T J Da d CT 

and, for all a G £ cx t> if f € £#;: 

M 2 ^ G m^hr I l^(")(^)| 2 dz, with e CT = (Vu) K • n K , CT - ZH^O ( 64 ) 

Since |i2jf| < — - ^ — |e CT | (where i?x is defined in (|25|0 . using the Cauchy-Schwarz 
Inequality, (|63|) and (|H4|) lead to the following bound: 

R 2 K < E m(o-)diir )CT ^ m(cr)dK, CT e^ < 

and, since dft- iCT < d CT and #£>>#: 

( < ^1 ^ / i^)wi 2 d- 

Then, (|27j) becomes: 

E KVpib^x - (Vu) K \ 2 m(K) < 



which gives the existence of C25, only depending on d and such that: 

E \(VvPdu)k ~ (Vu) K \ 2 m(K) < C^fa\\u\\ 2 H2[a) . (65) 

We have now to obtain an inequality similar to (|28|) (but without using u G C 2 (f2)). We will use 
here the fact that djf CT > 0diam(i<f) if a G 

If (j is a convex, bounded, open subset of M d , the well-known "Mean Poincare Inequality" gives, 
for all v G H 1 ^): 

\v(x)-muv\ 2 dx<-^dlm(B((),du)) [ \Vv(x)\ 2 dx, (66) 
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where m^v) is the mean value of v on uj, d w is the diameter of u, B(a,S) is the ball in Mr of 
center a and radius 5 and m(u>) (resp. m(B{a,5)) is the d-dimensional Lebesgue measure of uj 
(resp. B(a,5)). (A discrete counterpart of is given, for instance, in |Sj, Lemma 10.2.) 
Let K G M. We will use (jfifij) for cj = K. Since is the distance between xk to cr (for 
a G £r), there exists a G £k such that B(xk, d,K,a) C A. Then, one has 771(5(0,1))^ = 
m{B(xKi dx,a)) — m(if) and, using aV l( 7 > ^diam(A'), one obtains: 

m(A-) > 771(5(0, l))(0) d (diam(iT)) d . (67) 

Taking u = K in (|6"6"|). gives, for all K G ./VI and all v G H l (K): 

[ \v(x) - m w v\ 2 dx < -^diam(K) 2 / \Vv(x)\ 2 dx, (68) 
ift- # Jit 

Taking v equal to the derivatives of u (which are in H l (K) for all K G M) in (|68|) gives the 
existence of C26, only depending on d and 0, such that: 

Y, [ \Vu(x) - (Vu) K \ 2 dx < C^ 2 v \\u\\ 2 H2{n) . (69) 

Then, we conclude as in Lemma 1231 using (|B3|) and (|69j). that there exists Cj^] only depending 
on $7, # and u such that (|59|) holds. □ 

Theorem 4.2 (i/ 2 error estimate) Assume hypotheses 0^-^ d < 3, and that A and a 
are of class C 1 on 0. LetT) be an admissible finite volume discretization in the sense of Definition 
\2.1V and let 9 G (0, 6-p]. We assume that that card(£ft;) < 3 for all K G A4. Let u-p G Hx> be 
the solution of \34\j and u G Hq(Q) be the solution of J^J). We assume that u G H 2 (Q) (which is 
necessarily true if Q is convex). 

Let us first assume that Hypothesis \44\) holds. Then, there exists C27, only depending on Q, 9, 
a, (3, A and H^Hi^m), such that: 

\\uv - Pvu\\t> < Cffifiv, (70) 



\\uy - u\\ L 2 {n) < Cffifiy, (71) 

and 

\\V v uv - Vu\\ L2{n)d < Cffifiy- (72) 

(Recall that By, Vx> and \\ ■ \\<p are defined in Definition ^. SI Pp is defined in MU\) .) 

Let us then assume that \44\) no longer holds, then there exists C28, only depending on Q, 9, a, 

j3, A and \\u\\ H 2^, such that \70\ ), \71\) , \72\j hold with Gfegf/hy instead of Cffi^y- 

Proof. 

The proof of Theorem 14.21 follows the proof of Theorem 14.11 The quantities GfeS\ an d ( %IEi 
depending on 9, are now used to get a bound for Tfc^y) (as in [1]), and the quantity a ^ so 
depending on 9 since it is obtained with 1(59(1 (Lemma 14. 1(1 instead of ((21(1 (Lemma 14. 1(1 . is used 
to obtain a bound for 2^j(f). 

In order to obtain a bound for T^r^ v) (and then to conclude the proof of Theorem 14. 2|) . we need 
to obtain an inequality similar to (where the C 2 -regularity of u was used), which gives a 
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bound for the difference between the mean values of (A(x) — a(x)I^)Vu on K and on a if a G £k- 
Here, we will obtain a bound for the difference between these mean values using once again the 
consequence ()62|) of Inequality (9.38) in Lemma 9.4 of jHj. Applying ()62j) to the derivatives of 
(A — aId)Vu, there exists C29 only depending on $7, 8, A and a (indeed, the C 1 -norms of A and 
a), such that, for all K e M, all a e £ K and all v G H 1 ^): 

wk-m < q%j] m((j) (73) 

Following the proof of Theorem 14. 11 (|73|) is used to obtain a bound for r %^ y)' 

\^[j£v)\ < ^2 Yl \Vk\l - Vk\\A k ,l (v l -v k )\+ ^2 W - VkWAk^vkI J < 

|Mct - (J>k\ + \fJ>a ~ ML I / \ 1 <*cr« , IM(J-MaH / \ 7 $aV 

22 + £ m (-)^' 

where, in the last term, K is such that a G £k and where <5o-v = \vk — vl\ if <t = K\L G £ in t 
and = \vk\ i£ a =£ £ cxt ^£k- (We also used the fact that |Art,l| < and |Ak>| < ^jp-, 
thanks to 6t> > 0.) 

Then, using Cauchy-Schwarz Inequality and Q73|) . one obtains: 

< ^— ( ^ d CT (diam(A')||n||^ 2(x) + diam(L)||n||^ 2(L) ) 

o-=Jsr|Le£ int 

+ ^ d a di&m(K)\\u\\ 2 H2{K) )2 ■ 

0"G£ ex t 

Using ^ < 2hx>, diam(A') < hx> and the fact that card(£A') < \ for all K G M, one deduces 
the existence of Cq, only depending on ft, 8, A and a, such that: 

\^E$y)\ < C e h v \\u\\ H 2 {n) \\v\\ v . (74) 

Then, we conclude the proof of Theorem l4.2l exactlv as in the proof of Theorem l4.1l ( (|74|) replaces 
(PI)). □ 



5 Numerical results 

The scheme was tried for various academic problems, for which the analytical solution is known. 
For the Laplace equation, we compared the classical cell centered scheme to the new scheme, 
which we shall call the gradient scheme in the sequel. First note that in the classical cell 
centered scheme, the equation relative to a given cell involves the neighbors of this cell, while 
in the gradient scheme, it involves the neighbors of this cell and the neighbors of the neighbors. 
Hence in the case of a rectangular (resp. parallelipedic) mesh, the classical cell centered scheme 
is a 5 points (resp. 7 points) scheme, while the gradient scheme is a 13 points (resp. 24 points) 
scheme, Similarly, if one uses a triangular (resp. tetrahedral) mesh the classical scheme is a 4 
points (resp. 7 points) scheme, while the gradient scheme is a 10 points (resp. at most 17 points) 
scheme. Hence the gradient scheme is more expensive in terms of time and memory, although 



23 





Case 1 


Case 2 




homogeneous anisotropic 


heterogeneous anisotropic 




Rectangles FV 13 


Triangles VF10 


Rectangles FV 13 


Triangles VF10 


u 


2.00 


2.0 


2.2 


2.0 


Vu 


1.00 


1.0 


1.4 


1.3 



Table 1: Rates of convergence of FV13 and FV10 in a homogeneous anisotropic case and in a 
heterogeneous anisotropic case 



this is not so much, for example compared to the use of a Q finite element in the case of a 
parallelipedic mesh, which leads to a 27 points scheme. 

We tested the gradient scheme for some real anisotropic problems, the number of cells varying 
from 100 to 6400 in the rectangular meshes case (in fact, rectangles are squares), and from 700 
to 17500 in the triangular meshes case. The convergence rates have been computed by fitting a 
less-square regression on the logarithmic values of the errors and of the characteristic size of the 
mesh. 

The first case is an anisotropic homogenous problem with diffusion matrix 

( 1-5 0.5 \ 
A= U-5 1-5 J' 

The second case is a rotating permeability field, that is, the diffusion matrix is constant in the 
(r, 6) coordinates and equal to A r ^g = f ^ ^ J . The exact solution is taken to be u(x±, x 2 ) = 

\ ln((sci - .5) 2 + (x 2 - l.l) 2 ), on the domain Q =]0, l[x]0, 1[. The orders of convergence which 
were found are given Tabled 

Next, we tested different values of a to see how it affected the discretization error, on the first 
anisotropic case. Although the value of a does influence the resulting discretization error, the 
optimal value seems to be independent on the mesh, in both the triangular and rectangular cases, 
see Figure |2 Note that in the case of the error on the solution itself, the numerical optimal 
values for a are beyond the interval of convergence assumed in the theoretical analysis (0, 1). 
These numerical tests therefore indicate that this use of an discrete gradient in finite volume 
schemes leads to a correct numerical behavior, indeed comparable with low degree finite element 
schemes on similar problems. 

Finally, we replaced the point xk by the center of gravity of cell K in the definition (|16 )L I|17 |) 
of the coefficients Ak,l- I n this case, we recall (see Remark 12.3(1 that we obtain the discrete 
gradient based on the generalized Raviart-Thomas basis functions of ^0]. Indeed, the tests 
performed with this scheme for Case 1 or Case 2 did not yield correct approximations of the 
solution nor of its gradient. 

6 Conclusion 

In this paper, we constructed a discrete gradient for piecewise constant functions. This discrete 
gradient revealed several advantages: it is easy and cheap to compute, and it provides simple 
schemes for the approximation of anisotropic diffusion convection problems. We showed a weak 
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htb 




Figure 2: Diagrams of the errors on the solution (left) and its gradient (right) for various sizes 
of triangular (up) and rectangular (bottom) meshes, with respect to the value of the parameter 
a 
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property convergence of this discrete gradient to the gradient of the limit of the considered 
functions, together with a consistency property, both leading to the strong convergence of the 
discrete solution and of its discrete gradient in the case of a Dirichlet problem with full matrix 
diffusion. 

Since this notion of admissible mesh includes Vorono'i meshes, which are more and more used 
in practice, and which seem to remain tractable even in high space dimension, applications to 
financial mathematics problems are being studied [!]• Applications to finite volume schemes 
for compressible Navier-Stokes equations are also expected to be succesful j^H]- Further work 
includes a parametric study, and the generalization to meshes without the orthogonality condi- 
tion. 
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